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Background 


Why Quantifying Change in Bone via Bone Remodeling is Objective of NASA Digital Astronaut Project (DAP) 

• One of the main objeetives is to provide a tool to help HHC address Bone Gap Osteo 4 : We don’t know the eontribution of 
each risk factor on bone loss and recovery of bone strength and which factors are the best targets for countermeasure 
application; and Osteo? : We need to identify options for mitigation of early onset osteoporosis before, during, and after 
spaceflight. 

■ Skeletal loading along with endocrine regulation and local biochemical mediators are what drives the cellular 
mechanism of bone remodeling to maintain bone. 

■ Exercise induced loading, with appropriate input to a model can approximately predict the effect of specific exercise 
prescription and thus help to evaluate its benefits as a countermeasure option. Integrates with DAP Biomechanics 
Model and the DAP Muscle Model. 


Importance for the New Finite Element Based Strength Standard 

• Other main objectives intend to inform the HHC Bone Discipline's efforts to address Bone Gap Fracture 3 . We need a 
validated method to estimate the Risk of Fracture by evaluating the ratio of applied loads to bone fracture loads for 
expected mechanically-loaded activities during and after a mission 

■ One effort is underway to evaluate Finite Element (FE) estimates of bone strength (aka bone fracture loads) as a 
potential standard for bone health. 

■ A bone remodeling formulation that quantifies dynamic changes in bone has the potential of tracking changes in 
volume fractions that can relate to QCT BMD and ash density estimates, upon which FE bone strength is based [1]. In 
addition coupling a BR model with a QCT based FE model may also provide geometry changes. 


Bone Remodeling Model Implementation Plan 
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General Description of the DAP Bone Remodeling Model 


What does it do? 

It tracks changes in the bone when the balance between formation and 
resorption in the bone remodeling process becomes unbalanced. 

How does it do it? 

The cellular physiology, remodeling unit mechanisms, and mechano- 
transduction theory that drive the process are described mathematically. 

How does the computational algorithm work? 

Rates of change of bone volume fraction and cell populations are set to zero 
(Balanced healthy state with steady bone density). 

Balance is broken by skeletal unloading, and rate of change is no longer 0. 

The system including bone properties and cell populations are integrated in 
time to estimate the change. 



NOTE: Model parameters and methodology are currently focused on the femoral neck. 




Mathematical Description 


System of ordinary differential equations 


Bone \ oliune 
Fraction 

at Bq L q 

Eq. I 

Base Equation 

Osteoid 
\ oliune 
Fraction 

dO B C 

dt ^ B^ C„ 

{ ^ 
\0 + M^ 

\ 

-rO 

m 

/ 

Eq. 2 

Mineralized 

W^liime 

Fraction 

dM ^ C ( 

= r O — r — 

dt Q 1 

M ^ 

P + M J 


Eq. 3 

Responding 

Osteoblasts 

• Ptgf • ((1 Ptgf ) + PpGE ) • P,- 

Eq. 4 

Acthe 

Osteoblasts 

dR 

= D ((\ E ) + E ) B k -B-(\ E i 

VV^ ^TGF'^^PGE) ^ ^ PTH ^ 

at 

Eq. 5 

Acthe 

Osteoclasts 

=D -E A -E C 

^Cp ^RL ^TGF ^ 

Eq. 6 


State Variables and Definitions 


BVF{t) 

Bone Volume Fraction 

MO 

Cross Sectional Bone Area Resorhed per BRU 
(null-) 

M{t) 

Mineralized Volume Fraction 

Apit) 

Cross Sectional Bone Area Formed per BRU (imiF) 

0{t) 

Osteoid Volume Fraction 

1 

Actu ation DensiF in Normal State (#BRUs 
acti\ ated per da\' per nun2 ) 

Bit) 

Osteoblast Population (pM) 

A 

(Aj^ 'fa) Resorption Rate per Normalized 
( isteoclast Population 

C{t) 

Osteoclast Population (pM) 

h 

{Ap • ) Formation rate per Normalized 

Osteoblast Population 

Bo 

Reference (^isteohlast Population (pM) 

r 

m 

Mineralization Rate 

Co 

Reference Osteoclast Population (pM) 



Symbol Definitions in the Cell Equations 


B,. 

Concentration of Responding (')steohlasts 

Ac 

Rate of Elimination of C(t) (Apoptosis) 

B 

Concentration of Actn e Osteoblasts 

E 

J^TGF 

TGF-beta Receptor Occupancy Ratio 

C 

Concentration of Acth e (Osteoclasts 

E 

^ PGE 

Prostaglandin PGE-2 Receptor (Occupancy 

u 

Differentiation Rate of (Osteoblast Precursors 

E 

^ PTH 

Paratln roid Hormone Receptor Occitpancy 


Differentiation Rate of Responding (Osteoblasts 

E 

^NO 

Nitric (Oxide effect on RANKL 

u 

Differentiation Rate of Osteoclast Precursors 

E 

^RL 

RANKL Receptor Occupancy Ratio 

kp 

Rate of Elimination of Bit) 



Expressions for Osteroprotegerin 
(OPG) , RANKE and the ligand 
reeeptor eomplexes are derived 
via mass balanee equations. The 
eomplete detailed set of eellular 
dynamies is a eonsiderable 
modifieation of the work of 
Eemaire et al. [2] and Pivonka et 
al. [3] with the addition of 
effeetors related to skeletal 
loading. 


Modeling the Influence of Skeletal loading 


Some likely intermediaries that enable sensor cells to trigger effector cells are NO and PGE-2 [51. 

Released by Osteocvtes and Osteoblasts under mechanical stimulation 


r 

PGE2< 


Mediates differentiation of osteoblasts 
induced by TGF-^ 


V Stimulates proliferation of osteoblasts 


{ Stimulates production of OPG 
Inhibits production of RANKL 


The model gages the level of expression of NO and PGE 2 according to the level of bone apposition or 
bone resorption suggested by the daily strain s in Frost’s Mechanostat Theory as outlined below: 


Sensing strength or response level (SL) 
defined in relation to bone strain 




Complete Unloading £ = 0 SL = 0 
Remodeling Balance £ = £q SL = 1 



NOTE: Osteocytes are generally understood to be the sensor cells 


Production rate of NO and 
PGE 2 per cell are defined to 
be proportional to SL 


PGE Production Rate 
NO Production Rate 


Rate per cell 
Sq= PqXSLxY^x BVF 



= Pf^x SLx Fj X BVF 

t 

Rate per cell 


Osteocytes 

density 


Model Representation of Bone Loss Due to Insufficient Loading 






eiofw 
Physiological 
Maintenance 






occupancy ratio 
decreases 


RANKL increases OPG decreases 


j... J? 


imbaiancea m 

favor of 
resorption 







meraiized VO 
Fraction 
decrease 




Density (BMD) 
decreases 


Converting Experimental Data to Model Variables 


Definitions 


A Method for Mapping vBMD to BVF 


Ash density 


ash mass 
total vol 


Apparent (dry) density = 


dry bone tissue mass 
total vol 
ash mass 


Ash fraction a - — 

p^pp inorganic mass + organic mass 

- density Mineralizd bone - density of Osteoid 


Given: 

A Pre Bed Rest QCT BMD value. 

A Bed Rest Duration Length of N days. 

A Post Bed Rest QCT value. 

1. (a) Convert Pqct^^ Pash (^-S- K^y^k regression) 

(b) Convert p^^^ to p^^^ (e.g. Schileo regression) 

(c) Compute initial ash fraction a = p^^^ / p^^^ 


Ash fraction has a theoretical limit 0.7 [4]. 
apparent density 


BVF = 


true tissue density 


p^{glcc) = \Al + \.29a [4] 


a = 


MD^ + OD^ 


M D„ -\-0'D 


or 


m 


o 


M + O 


= p^ (Used in Model) 


2. Initial value M = / (0.7 x D„) 

Solve for initial value O using a definition. 

3. Run computational simulation subject to loading 
history (i.e. bed rest) for N days to track change 

in M, O, a, p^^^, Pqqt^ (BMD), md BVF 

4. Compare BMD to QCT BMD 
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Preliminary Validation Results for Bone Deconditioning Simulations 
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Preliminary validation results for predicting: 

(a) Group mean BVF [9-11]. 

(b) Loss of trabecular bone after 70 days of bed rest. 

(c) Loss of cortical bone after 70 days of bed rest, 

(d) Time course change of mean DXA BMD for 18 control subjects during 17 weeks of bed rest [12]. 


The QCT bone analysis data was provided by the NASA Johnson Space Center Bone Lab through the NASA Life Sciences Data Archive 


Future Work 


Near Term: 

• Develop/formulate a daily load formula for quantifying exercise 
induced loading and test against exercise treated subjects (e.g. 
CFT70 study) 

Long Term : 

• Develop method for transforming force data from biomechanics 
modeling of specific exercise devices into stress/strain input 

• Integrate the computational model with Finite Element Method 

• Validate model using QCT data from spaceflight research 

• Develop model for predicting bone adaptation for trochanter, total 
proximal femur and lower lumbar 

• Bone adaptation prediction for more than 1 80 days of spaceflight 
exposure with exercise countermeasure 
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